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Preface 


The  ability  to  control  the  ringing  effects  of 
Fresnel  diffraction  is  a  necessity  in  high  power  laser  sys¬ 
tems.  Apodisation  is  a  technique  employed  in  this  study  to 
combat  this  phenomenon. 

Apodisation  is  a  topic  of  current  interest  in  the 
optics  community  and  certainly  to  the  Air  Force  with  the 
applications  possibility  being  employed  to  weapon  systems. 

I  wish  to  thank  Major  James  P.  Mills  for  proposing 
the  topic  and  giving  such  helpful  advice  and  direction; 
also,  Mrs.  Pat  Mills  for  her  support  in  the  production  of 
the  many  computer  plots.  Additionally,  I  would  like  to 
thank  Mrs.  Phyllis  Reynolds  who  did  such  an  excellent  job 
in  typing  this  thesis;  as  well  as  Mr.  Ron  Gabriel  for  his 
outstanding  support  in  the  laboratory.  Above  all,  I  would 
like  to  thank  my  family;  my  children  Tom  and  Erica,  and  my 
wife  Terese  for  their  love,  patience,  and  understanding. 
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Abstract 

The  effects  of  Gaussian  apodisation  with  -0.8  waves 
of  spherical  aberration  on  a  coherent  optical  system  was 
examined.  The  area  of  study  was  in  the  near-field  concen¬ 
trating  on  the  areas  near  the  focus  in  the  Fresnel  region 
of  an  f/12  system  using  optical  wavelengths.  Computer  pre¬ 
dictions  were  made  for  four  cases:  unapodised-unaberrated, 
apodised-unaberrated ,  unapodisec-aberratea ,  and  apodisea- 
aberrated.  Predictions  were  made  for  cross-sectional  planes 
perpendicular  to  the  optical  axis  using  Fourier  optics. 
Meridional  plane  predictions  were  produced  using  a  numerical 
integration  method  of  determining  the  Kirchoff  integral. 
Additionally,  experimental  data  are  given  to  compare  with 
the  predictions.  It  is  shown  that  the  experimental  data 
matches  the  computer  predictions  and  that  apodisation  is  an 
effective  method  for  controlling  the  ringing  due  to  edge 
effects  and  spherical  aberrations.  Additionally,  fine 
structure  corresponding  to  Young's  double  slit  interference 


THE  NEAR-FIELD  EFFECTS  OF  APODISATION 


ON  COHERENT  ABERRATED  OPTICAL  SYSTEMS 


I.  Introduction 


This  thesis  examines  the  effects  of  apodisation  in 
the  near-field  of  a  coherent  optical  system  with  aberra¬ 
tions  both  experimentally  and  theoretically.  Apodisation 
is  the  intentional  alteration  of  the  amplitude  transmittance 
of  an  optical  system.  It  is  known  that  apodisation  will 
effectively  smooth  the  near-field  irradiance  fluctuations 
caused  by  diffraction  from  the  edges  of  apertures  in 
unaberrated  optical  systems.  Thompson  and  Krisl  (25:109) 
provide  graphic  evidence  of  this.  This  thesis  will  study 
the  near-field  effects  of  apodisation  on  aberrated  optical 
systems . 

The  method  undertaken  involves  theoretical  predic¬ 
tions  of  the  apodised  response  to  an  aberrated  system 
using  scalar  diffraction  theory.  This  is  done  by  computer 
simulations  which  model  the  near-field  irradiance  and  phase 
under  four  circumstances: 

1)  the  unaberrated  and  unapodised  near-field, 

2)  the  aberrated  near-field, 

3)  the  apodised  near-field,  and 

4)  the  aberrated  near-field  under  apodisation. 


In  this  manner  the  effects  of  aberrations  and  apodisation 
may  be  clearly  seen. 

These  predictions  are  then  compared  with  experi¬ 
mental  evidence  of  the  near-field  under  the  same  circum¬ 
stances  for  spherical  aberration. 

With  the  advent  of  high  power  lasers,  diffraction 
effects  can  cause  serious  damage  to  components  in  the  opti¬ 
cal  chain.  Control  of  these  effects  is  necessary  for 
efficient  and  safe  operation  of  high-power  optical  systems. 
This  study  takes  the  elementary  case  of  prim, ary  aberrations 
on  a  plane  wave  incident  upon  a  lens  and  aperture.  The 
resulting  diffraction  effects  on  the  converging  wavefront 
are  then  analyzed  in  the  Fresnel  region  with  and  without 
aberrations,  both  under  the  effects  of  apodisation  and 
without  apodisation. 

Figure  1  shows  the  effects  of  apodisation  on  an 
unaberrated  wave,  and  a  wave  with  third  order  spherical 
aberration  of  -0.8  waves.  Figure  1  shows  a  cross- 
sectional  graph  at  the  first  on-axis  minimum  in  irradiance 
from  the  focus  toward  the  aperture  in  an  f/12  system  with 
a  Fresnel  number  of  439.021  and  a  632.8  nm  wavelength. 

These  parameter  values  are  used  throughout  the  study  and 
arise  from  the  experimental  system  described  in  Section  IV. 
The  irradiance  is  normalized  to  one  at  the  focus  for  the 
unaberrated  unapodised  system.  The  radial  coordinate  is 
also  normalized  (see  Section  III). 


Background 


The  near-field  has  been  examined  by  several  authors, 
most  notably  are  Campbell  and  DeShazer  (4)  who  evaluated 
the  near-field  amplitudes  and  irradiances  using  the 
cylindrical  coordinate  representation  of  the  Fresnel 
approximation  to  the  Kirchoff  diffraction  integral.  This 
was  done  for  Gaussian  apodisation.  Olaofe  (22)  has  produced 
an  analytical  solution  to  the  same  integral.  Neither  of 
these  works  consider  aberrations. 

Harvey  and  Shack  (12)  have  studied  aberrations 
through  the  evaluation  of  the  Rayleigh-Sommerf ield  diffrac¬ 
tion  formula  generalizing  into  a  two-dimensional  Fou  ier 
transform.  This  method  yields  results  which  do  not  suffer 
from  the  paraxial  and  Fresnel  approximations.  This  work, 
however,  is  mostly  concerned  with  the  area  of  the  focal 
plane . 

Avizonis,  et  al.  (1)  studied  the  near-field  aber¬ 
rations  with  the  purpose  of  gaining  insight  into  their 
effects  at  the  focus.  The  evaluation  method  used  was  also 
by  the  cylindrical  Fresnel-Kirchof f  integral. 

The  study  of  phase  has  been  done  principally  by 
Farnell  (7;  8;  9)  who  did  his  research  using  microwaves  and 
the  calculation  methods  of  Linfoot  and  Wolf  (17)  for  theo¬ 
retical  evaluation.  This  was  done  for  both  unaberrated  and 
aberrated  cases. 
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Bachynski  and  Bekefi  (2)  examined  both  irradiance 
and  phase  on  microwave  systems  and  performed  a  detailed 
analysis  on  the  effects  of  aberrations.  Most  of  their  work 
considered  the  areas  at  or  near  the  focus.  Their  work  is 
also  quite  notable  in  their  prediction  and  measurement  of 
an  axial  irradiance  in  the  near-field  greater  than  that 
at  the  focus  for  microwaves. 

Work  specifically  on  apodisation  has  been  done  by 
Thompson  and  Krisl  (25)  and  Kuzimina,  et  al .  (15).  Both 
papers  evaluated  apodisation  effects  of  various  types  of 
apodisers.  The  latter  paper  is  concerned  much  more  with  the 
focal  plane.  The  former  paper  demonstrates  the  ringing 
decrease  in  image  formation  and  also  in  the  Fresnel  region. 

Holmes,  et  al.  (13)  examine  Gaussian  focused  beams 
using  numerical  integration  in  cylindrical  coordinates  of 
the  approximated  Kirchoff  equation. 

The  study  of  apodisation  on  near-field  aberrations 
is  a  significant  void  in  the  literature,  with  some  work  on 
near-field  aberrations  and  some  on  apodisation,  but  none  on 
both.  The  most  complete  study  of  apodisation  on  aberra¬ 
tions  in  the  image  plane  is  by  Mills  (21) .  This  thesis  is 
actually  an  extension  to  the  work  done  by  Mills,  in  that 
the  apodisation  effects  on  aberrations  are  now  being  taken 
into  the  near-field.  Additionally,  this  thesis  is  pri¬ 
marily  concerned  with  spherical  aberration,  both  theoreti¬ 
cally  and  experimentally.  Also,  this  study  considers  a 
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very  large  aperture  with  respect  to  the  wavelength  of 
light.  The  aperture  diameter  is  over  21,000  wavelengths 
wide. 


II .  Theory 


Near-field  Fourier  Optics 

The  field  amplitude  in  the  near-field  is  given  by 
solutions  to  the  Helmholtz  equation.  Scalar  diffraction 
theory  is  based  on  the  Green's  function  solution  to  this 
equation. 

An  integral  solution  to  the  amplitude  U(x^,y^,z) 
in  the  near-field  of  an  aperture  illuminated  by  a  plane 
wave  is  given  by  the  Kirchoff  diffraction  integral: 

U(x1>y;L,z)  =  fx//vU(X0,y0)  ~"s —  (1  +  cosX>  d°  (1) 

where  A  is  the  wavelength,  k  is  the  wave  number  2tt / A  and, 
U(x0,y0)  is  the  complex  field  amplitude  in  the  aperture. 

The  integration  is  over  the  wave  front  in  the  aperture 
area  l  using  the  geometry  in  Figure  2.  The  parameter  s  in 
Equation  (1)  is  the  line  segment  connecting  an  arbitrary 
aperture  point  Q  to  an  observation  point  P.  The  angle  X 
is  between  s  and  the  normal  to  the  incident  wave  front  n. 
For  an  unaberrated  system  the  incident  wave  is  a  plane 
wave  so  n  is  perpendicular  to  the  Xq  and  y^  axes  and 
parallel  to  z. 


Several  approximations  are  now  made  to  simplify 
Equation  (1) .  The  first  is  the  paraxial  approximation 


and  assumes  that  the  angle  x  is  small  so 


1  +  cos  x  =  2 


(2) 


Secondly,  s  is  approximated  to  be  nearly  equal  to  z.  There¬ 
fore,  s  is  replaced  by  z  in  the  denominator  of  Equation  (1) . 
This  will  not  suffice  in  the  argument  of  the  exponent  because 
the  error  in  the  approximation  would  be  greatly  magnified 
since  it  is  multiplied  by  k,  a  very  large  number.  The 
Fresnel  approximation  is  therefore  applied  to  s  in  the 
exponent.  The  exact  value  of  s  is  given  by 


s 


1  + 


2 


+ 


(3) 


Using  the  binomial  expansion  on  the  square  root;  keeping 
the  first  two  terms,  and  then  substituting  back  into  the 
Kirchoff  integral  gives 

oc 

u(xl'yl'z)  =  TTzJJ  U(W  exp  | zl  [(x0‘xl>2 

+  (yo_yi)2]|dxodyo  (4) 

The  limitations  due  to  the  approximations  are  the 

most  severe  in  the  region  closest  to  the  lens.  The  paraxial 
approximation  has  an  accuracy  good  to  within  5  percent  for 
an  angle  x  up  to  18°.  The  Fresnel  approximation  is  accurate 
to  the  limits  of  the  paraxial  approximation.  The  approxi¬ 
mations  limit  the  closest  approach  to  the  aperture  to 


12.7  percent  of  the  focal  length.  This  is  as  measured  on 
the  optical  axis  from  the  aperture  for  an  f/12  system. 

Equation  (4)  is  recognized  as  a  convolution  integral 
and  is  quite  useful  in  many  applications  of  diffraction 
theory.  For  this  thesis  it  is  most  useful  in  a  form  found 
after  rearranging  the  exponentials: 


U  (x^y^z) 


exp 

iXz 


i2  ttz 
X 


exp 


in  (x2 

Xz  'X1 


+  Yi)]i/1u  <X( 


'V 


*  exp[^  lxo +  yo>]exp[icr  ,xoxi +  yoyi']|dxodyo  ,5) 

where  U(Xg,yg)  defines  the  complex  field  amplitude  in  the 
aperture  and  is  given  by 


U(xO'yo’  =  Uin<X0'y0)T<X0-y0,G<X0'y0)A<X0'y0)L<x0'y0)  (6) 

where  U^n  is  the  incident  amplitude,  T  is  the  aperture 
transmittance,  G  models  the  apodisation,  A  represents  the 
aberrations  occurring  in  the  aperture,  and  L  is  the  dis¬ 
turbance  imparted  on  the  field  due  to  a  lens.  Each  of 
these  terms  will  be  looked  into  in  detail  later. 

Equation  (5)  is  a  two-dimensional  Fourier  transform 
of  the  product  of  the  aperture  function  U(xQ,y0)  and  the 
quadratic  exponential  term  which  describes  the  near-field 
region  dependence  of  the  integrand.  This  latter  term  will 
be  called  the  Fresnel  term.  Another  form  for  Equation  (6) 
is 
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U(x1>Yl,z)  =  ^  exp[i^5]exp 


x  F  U(x0,y0)  exp  il  (x*  +  y*) 


[ff  (x?  +  yl'] 
l  lx0  +  y0»]| 


where  F { }  indicates  the  Fourier  transform  operator.  The 


transform  spatial  frequencies  are  f  =  x /xz  and  f  =  y  /xz 

xu  y  u 


The  irradiance  is  given  by 


Kx^y^z)  =  u*(x1  fyn  FZjutx,  ,y,  ,z) 


1,J1 


1  J 1 1 


where  the  asterisk  denotes  the  complex  conjugate.  This 


equates  to 


I(x1,y1,z)  =  — i-j  Fju(xo,y0)  exppglxj  +  yj 


The  modulus  is  then  simply  the  square  root  of  the  irradi¬ 


ance  : 


M(xi,yi,z)  =  T£  Fju(xo,yo}  exp[S(x 


2^2. 

0  +  V 


Cylindrical  Coordinate  Representation 


If  the  aperture  function  U  (xQ,y0)  is  circularly 


symmetric  then  the  following  substitutions  apply: 


,  2  ^  2, 
p  =  <x0  +  y0> 


Xq  =  p  cos  e 


0  =  arctan (y0/xQ) 


yu  =  P  Sin  e 


2  2 

r  =  (f  +  f  ) 
x  y 


f  =  r  cos  f 
x  T 


i>  =  arc  tan  ( f  /  f  ) 
y  * 


f  =  r  sin  4i  (11) 


mmism 


Substitutions  are  made  into  Equation  (5)  and  the 

following  trigonometric  Bessel  function  identity  applied 

2  TT 


J0  (b)  =  ~  /exp  [-ibcos  (6  -41)]  de 
71  J0 


(12) 


where  b  is  any  substitution  variable.  The  irradiance,  in 
cylindrical  coordinates  is  then  given  by 


I  ( 


r ,  z)  =  J U(p)  exp(iiTp2/Xz)  JQ  (2frp/Az)  pdp 


(13) 


Equation  (13)  is  the  primary  equation  used  to  cal¬ 
culate  the  values  of  the  irradiance  both  on  the  optical 
axis  and  in  the  radial  direction. 


On-axis  Solutions .  On  the  optical  axis,  Equation 
(13)  simplifies  since  r  goes  to  zero  and  JQ{0)  =  1.  U(p) 
may  still  be  quite  complicated  so  this  simplification  does 
not  always  lend  itself  to  analytical  solutions.  It  has 
been  solved  for  the  case  of  an  unaberrated  unapodised  wave 
and  can  be  integrated  directly  yielding 


I  (0  ,  z ) 


(f1 

2 

sin  (<pQ/2) 

,zl 

Uq/21  ! 

(14) 


where  4>g  is  the  term  arising  from  the  combination  of  the 
Fresnel  term  and  the  term  due  to  a  lens  of  focal  length  f 
placed  in  the  aperture: 


t 


0 


(15) 
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where  A  is  the  area  of  the  aperture.  Olaofe  (22:1654)  and 
Mahajan  (19:8)  have  both  solved  for  the  case  of  Gaussian 
wave  fronts. 

For  the  cases  of  spherical  aberration  and  apodisa- 
tion  with  spherical  aberration,  the  integral  becomes 
analytically  unmanageable  and  requires  numerical  methods 
for  solution. 


The  Aperture  Function 

The  Huygens-Fresnel  principal  states  that  the  field 
amplitude  at  any  point  is  given  by  the  superposition  of  an 
infinite  number  of  spherical  waves  radiating  from  an 
infinite  number  of  points  in  the  aperture.  Thus,  the  aper¬ 
ture  function  describes  the  scalar  field  amplitude  in  the 
aperture.  This  is  given  by  Equation  (6): 


U(x0'y0)  = 


U.  (x. 
m  0 


-yo)T|xo 


,y0)G(xQ,yc 


)A(x0'y0)L(xn'yn) 


0,y0‘ 


(16) 


Each  term  will  be  examined  one  at  a  time. 


Incident  Wave  (U .  ).  This  term  is  simply  the 
-  -  — in— 

complex  field  amplitude  incident  upon  the  aperture.  In 
this  study  this  will  always  be  considered  as  a  plane  wave 
of  unit  amplitude. 

Transmittance  Function  (T) .  The  transmittance  func¬ 
tion  defines  how  the  boundary  conditions  of  the  aperture 
allow  the  incident  infinite  plane  wave  to  transmit  through 
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the  aperture.  For  the  case  of  a  circular  aperture  this  is 
a  circle  function  of  unit  amplitude  out  to  the  aperture 
boundary: 

!1  0  _<  p  £  a 

(17) 

0  otherwise 

Apodisation  (G) .  With  a  Gaussian  apodiser  the 
transmittance  is  a  Gaussian  function  truncated  by  the  edges 
of  the  circle  function.  The  apodiser  could  be  a  complex 
function  but  this  study  will  consider  only  real  Gaussian 
apodisers  for  three  reasons: 

1)  the  Fourier  transform  of  a  Gaussian  is  itself  a 
Gaussian, 

2)  the  impulse  response  is  real  and  positive,  and 

3)  there  are  strong  similarities  with  some  laser 
wave  fronts . 

The  equation  for  a  Gaussian  is 

G  ( p)  =  exp(-2yp2)  (18) 

2  -2 
where  y  =  (a/w)  ,  w  being  the  irradiance  e  beam  radius. 

In  this  study  y  =  3/2.  Often  the  coefficient  2>  is  leplaced 

by  1/^q  in  the  literature.  In  either  case  the  apodisation 

function  considered  here  has  the  form 

G ( p )  =  exp(-3p2)  (19) 


This  function  has  the  property  of  having  a  transmittance  of 
0.05  at  the  aperture  boundary  thus  only  weakly  truncating 
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the  wave  front  as  in  Figure  3.  The  weak  truncation 
eliminates  almost  all  the  edge  effects  as  seen  in  the  work 
by  Mills  (21:3). 


Aberrations  (A) .  The  aberration  of  primary  impor¬ 
tance  in  this  study  is  third  order  spherical  aberration. 

It  is  represented  mathematically  by  the  Zernike  polynomials 
and  monomials  which  represent  the  aberrated  wave  in  the 
aperture.  A  brief  discussion  of  the  Zernike  polynomials 
follows  which  parallels  the  treatment  by  Malacara  (20:489- 
505)  and  Born  and  Wolf  (3:464-466)  which  should  be  con¬ 
sulted  for  greater  detail. 

The  Zernike  polynomials  are  orthonormal  terms  which 
are  balanced  polynomials  that  provide  the  highest  Strehl 
ratio  at  the  paraxial  focus.  The  poxynomials  satisfy  the 
condition 

1  2v 

I  I  z£*z%  dp  dG  =  — jy  c  (20) 

JJ  n  m  n+1  nm 

0  0 

where  6  is  the  Kronecker  delta  function  and  the  asterisk 
nm 

£ 

represents  the  complex  conjugate.  The  function  z^  can  be 
broken  into  a  radial  term  and  an  axial  term  or 

Z*  =  Rl  (p)  exp  (  i £  6 )  (21) 

n  n 

where  n  is  the  degree  of  the  polynomial  and  £  is  the 
angular  dependence.  The  radial  distance  p  is  normalized  to 
the  aperture  radius.  Both  n  and  >.  must  either  be  odd  or 
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0 

even.  A  radial  polynomial  Rn  exists  for  each  pair  of 
numbers  n  and  | £ | .  Also 

r£  =  r'£  =  rI  1 1  (22) 

n  n  n 

When  n  is  even  the  polynomial  is  symmetrical  and  the 
exponents  of  the  polynomial  terms  are  even.  When  n  is  odd 
the  polynomial  is  antisymmetric  and  the  exponents  are  odd. 
The  radial  polynomials  are  found  by  the  relation 


m 


R^'2m  (p) 


(-1) 


(n-s)  1 


n-2  s 


s!  (m-s)  !  (n-m-s) 


s=0 


(23) 


For  defocus,  n  =  2  and  £  =  0.  For  third  order  spherical 
aberration,  n  =  4  and  £  =  2.  Since  n  is  even  in  both  cases 
there  is  no  angular  dependence. 

The  Zernike  polynomials  may  be  converted  to  rec¬ 
tangular  coordinate  monomial  forms  by  the  substitutions  of 
x  =  psinG  and  y  =  pcosG .  The  Zernike  polynomials  and 
monomials  are  displayed  in  Table  I  which  is  reprinted  from 
Malacara.  Only  the  terms  which  correspond  to  primary 
aberrations  have  been  kept.  Coma  and  spherical  aberrations 
listed  are  third  order. 

When  the  polynomials  are  included  in  the  aperture 
function  they  are  in  the  argument  of  a  complex  exponential 
to  represent  the  phase  disturbance  on  the  wavefront.  They 
are  also  multiplied  by  the  amplitude  of  the  aberration 
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TABLE  I 


ZERNIKE 

!  POLYNOMIALS 

AND  MONOMIALS 

(20:493) 

n 

m 

n-2m 

Polynomials 

Monomials 

Meaning 

0 

0 

0 

1 

1 

Constant  term 

1 

0 

1 

psinO 

X 

Tilt  in  x 

1 

-1 

pcos8 

y 

Tilt  in  y 

2 

0 

2 

2 

p  sin2  6 

2xy 

Astigma.  ±45° 

1 

0 

(2p2-l) 

2  2 
-l+2yz+2x 

Focus  shift 

£. 

-2 

2 

p  cos2  9 

2  2 
y  -x 

Astigma.  0°,90° 

3 

1 

1 

(3p3-2p) sin6 

-2x+3xy2+3x3 

x  axis  coma 

2 

-1 

(3p3-2p) cos9 

-2y+3y3+3x2y 

y  axis  coma 

4 

2 

0 

6p4-6p2+l 

2  2  4 

l-6y  -6x  +6y 

2  2  4 

+12x  yZ+6x^ 

Spherical 

in  wavelengths.  Plots  of  the  Zernike  polynomials  are  in 
Figure  4. 


Lens  Term  (L) .  The  inclusion  of  a  converging  lens 
creates  a  spherical  wavefront  with  its  center  at  the  focus. 
The  phase  disturbance  for  the  lens  is  represented  by 

L(p)  =  exp  (— in p2 / Xf )  (24) 

Equation  (24)  is  a  quadratic  approximation  to  the  actual 
spherical  wavefront.  This  representation  for  the  lens 
effects  is  pervasive  in  the  literature  and  is  developed  by 
Goodman  (10:80)  as  well  as  other  sources.  A  converging 


lens  is  always  considered  to  be  in  the  aperture  for  all 
cases  in  this  thesis. 
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III.  Computer  Model inc 


Programs 

The  prediction  of  the  irradiance,  modulus,  and 
phase  at  various  points  in  the  image  field  of  an  aperture 
was  done  by  computer  modeling.  Three  main  approaches  were 
used  and  are  described  below. 


Fourier  Transform  Method.  Values  for  the  irradi¬ 
ance,  modulus,  and  phase  can  be  obtained  for  cross- 
sectional  planes  by  performing  the  two-dimensional  Fourier 
transform  of  Equations  (9)  and  (10) .  This  code  is  based  on 
creating  a  two-dimensional  array  of  the  aperture  function, 
multiplying  by  the  Fresnel  term,  and  then  taking  the  Fourier 
transform.  The  code  was  initially  developed  by  Mills  (21) 
to  give  results  in  the  focal  plane.  Its  modification  now 
allows  computation  for  results  in  the  near-field.  The  code 
is  named  N256.FOR  and  is  included  in  the  appendix.  The 
Fourier  transform  is  calculated  by  the  IMSL  (14)  routine  FFT3D. 

To  be  generic  for  any  system.  Equation  (9)  was  nor¬ 
malized  by  the  variable  substitutions 


xn  '  x0/a 


x  =  X,  / (Az/a) 
n  I 


=  Va 


yn  =  Yl/  (Az/a) 


In  =  I0  <Vy0'z)/ (7ia/Af) 


m 


■S' 


This  yields  the  normalized  equation  after  dropping  sub¬ 
scripts 

I  / 

( 1— c ) 


K 


l  f(  i 

x^y^z)  =  - - jZ<t(x0,y0)exp  ittFN- 


,  2  ^  2. 

x  (xQ  +  yc) 


]) 


(26) 


where  c  is  simply  the  ratio  z/f  and  FN  is  the  Fresnel  num- 

2 

ber  of  the  system,  a  /Af.  in  this  study  a  Fresnel  number 
of  439.021  was  used  which  uses  the  parameters  of  the  experi¬ 
mental  system  with  A  =  6.328  x  10  5  cm,  an  aperture  radius 
of  0.66675  cm,  and  a  focal  length  of  16.002  cm.  The  reasons 
for  these  values  are  defined  by  the  experiment.  With  the 
paraxial  approximation  on  this  system,  accuracy  to  within 
5  percent  occurs  when  z/f  is  less  than  or  equal  to  0.12. 

The  areas  under  study  here  are  well  within  the  approxi¬ 
mations  used. 

Cylindrical  Coordinate  Numerical  Integration .  The 
near-field  irradiance  can  also  be  given  by  Equation  (13) 
which  uses  cylindrical  coordinates  and  the  zero  order 
Bessel  function  of  the  first  kind.  As  also  shown  by 
Mahajan  (18:3036),  it  is  normalized  with  the  substitutions 

Pn  =  P/a,  rn  =  r/(Af/2a),  IR(rn)  =  IQ/(^a2/Af)2  (27) 

Equation  (13)  now  becomes,  after  dropping  subscripts 

2 


I( 


r,z)  =  4^|  j  J  £exp  iF4  (o)J 


<*<> 


(28) 
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where  I<t(p)  is  the  sum  of  the  Fresnel  and  aberration  terms. 

Radial  data  can  be  computed  using  this  equation  as 
well  as  data  on  the  optical  axis.  The  codes  use  Simpson's 
rule  integration  and  the  Bessel  function  is  computed  by  the 
IMSL  routine  MMBSJO  (14) .  An  alternative  to  using  the  IMSL 
routine  is  to  use  the  Bessel  identity  (5:188) 

TT 

Jq  (b)  =  —  I  2cos  (bsin<j>)  dp  (29) 

v  J0 

using  numerical  integration  in  the  subroutine  to  calculate 
JQ  (b) .  The  values  in  the  meridional  plane  are  calculated 
by  this  method  by  the  code  MERID.FOR  which  is  also  included 
in  the  appendix. 


Double  Integration  Method.  The  double  integration 
method  follows  directly  the  method  used  by  Bachynski  and 
Bekefi  (2:430)  in  the  generation  of  their  data,  as  well 
as  the  works  done  by  Lee  and  Farnell  (16:273) 


-ika  U. 

U(r,z)  =  ■■ 2vzf  in  exp[ik(f 


0  0 


2v 

2 

exp (-Gp  ) 


x  exp[-ipp  -  iqpcos (6-v)  +  ikV(p, 6) ] pdpde 


(30) 


where 


p  =  (ka  /2zf) (f-z) 
q  =  (ka/ f ) r 
V  =  aberration  terms 
G  =  Gaussian  apodisation  constant 


Equation  (30)  utilizes  the  cylindrical  coordinate  system 
previously  mentioned,  but  does  not  utilize  the  Bessel 
identity  in  Equation  (12)  keeping  the  angular  dependence  in 
the  integral.  This  allows  the  computation  of  non-circular ly 
symmetric  aperture  functions  which  can  be  used  to  study  the 
effects  of  coma,  astigmatism,  and  tilt  aberrations.  This 
program  utilizes  the  IMSL  integration  routine  DBLIN  (14). 

This  method  provides  very  accurate  results  due  to 
the  cautious  Romberg  integration  method  used  by  the  IMSL 
routine.  A  great  drawback  to  this  method  is  the  enormous 
amount  of  computer  time  which  the  code  takes.  The  program 
which  utilizes  this  method,  IRAD.FOR,  is  also  included  in 
the  appendix. 

Computer  Predictions 

Studied  Regions .  There  were  two  main  regions  of 
the  near-field  which  were  closely  studied:  the  area  near 
the  focus  and  the  area  near  z/f  =  0.95,  as  shown  in 
Figure  5.  The  focal  area  was  studied  because  there  was 
existing  data  to  check  the  accuracy  of  the  codes'  predic¬ 
tions.  Also,  it  is  an  area  of  interest  in  that  it  provides 
a  unique  opportunity  to  see  the  physical  responses  in  the 
transition  from  Fresnel  to  Fraunhofer  diffraction.  All 
plots  were  made  with  the  graphics  software  package  DISSPLA 
(6)  . 
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The  area  near  z/f  =  0.95  was  studied  to  observe 
the  effects  of  aberrations  and  apodisation  in  the  Fresnel 
region.  Areas  closer  to  the  aperture  were  found  to  require 
enormous  computer  resources  far  beyond  the  capabilities  of 
the  systems  available.  The  creation  of  enormous  two- 
dimensional  arrays  is  required  to  accurately  predict  the 
cross-sectional  disturbance  in  these  areas.  Even  at  z/f 
=  0.95,  cross-sectional  arrays  were  incalculable  due  to  the 
required  size  for  accurate  measurement.  A  look  at  the 
photographs  in  the  next  section  demonstrates  the  vast  amount 
of  resources  required  to  predict,  store,  and  provide  an 
accurate  output  of  such  large  amounts  of  information. 
Fortunately,  radial  plots  in  this  region  could  be  produced 
without  too  much  difficulty.  Radial  plots  are  sufficient 
because  the  cross-sectional  plane  is  circularly  symmetric. 

The  Focal  Region.  The  focal  region  under  study 
consists  of  the  area  from  z/f  =  0.98  to  z/f  =  1.02  and  from 
rn  =  -4.0  to  4.0  using  the  normalized  radial  coordinates 

rn  =  r (2a/>  f )  (31) 

The  system  studied  throretically ,  used  the  param¬ 
eters  of  the  experimental  system  with  a  Fresnel  number  of 
439.021  . 

The  meridional  plane  over  the  focal  region  is  given 
in  Figure  6.  The  meridional  plane  is  the  plane  defined 
when  the  angle  ••  is  zero  as  seen  in  Figure  2.  It  can  be 
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Meridional  Plane  Irradiance  Near  the  Foe 
a)  Topographical;  b)  Contour  Plot  from 
Linfoot  and  Wolf  (17:826) 


easily  compared  with  the  two-dimensional  contour  plot  from 
Linfoot  and  Wolf  (17:826)  which  has  a  different  coordinate 
system  covering  a  slightly  smaller  area  along  the  optical 
axis,  but  the  same  area  on  either  side  of  the  axis.  This 
qualitatively  demonstrates  the  accuracy  of  the  code 
MERID.FOR  which  produced  the  plotted  data. 

The  effects  of  apodisation  and  spherical  aberration 
on  the  entire  focal  area  in  the  meridional  plane  can  be 
seen  quite  dramatically  in  Figure  7.  The  spherically 
aberrated  modulus  changes  the  values  seen  in  the  unaberrated 
plot.  Phase  is  plotted  in  Figure  8.  The  peak  in  modulus 
corresponds  to  the  circle  of  least  confusion.  When  the 
apodiser  is  in  place,  the  ringing  in  the  unaberrated  system 
completely  disappears  in  this  region.  In  the  aberrated 
system,  the  ringing  is  very  much  reduced  but  not  completely 
eliminated.  In  all  cases  of  apodisation,  there  is  a  sig¬ 
nificant  reduction  in  amplitude  of  the  peak  values.  This 
is  due  to  the  attenuation  of  the  energy  allowed  to  pass 
through  the  aperture  by  the  absorbing  apodiser.  Addi¬ 
tionally,  the  edge  effects  have  been  dramatically  reduced, 
since  the  edges  are  only  5  percent  as  high  as  in  the 
unapodised  aperture. 

Many  of  the  two  pi  jumps  in  phase  are  due  to  a 
recycling  of  the  arctangent  function  used  to  calculate  the 
phase  value.  A  false  impression  of  a  phase  jump  is  thus 
given.  A  graphic  indication  of  the  rapid  phase  changes 
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Meridional  Plane  Phase  in  Radians 
Near  the  Focus 


occurring  in  the  aberrated  plots  prior  to  the  focal  plane 
is  easily  seen  however.  The  phase  shifts  behind  the  focal 
plane  are  real  and  can  be  seen  to  represent  the  diverging 
rings  in  the  modulus  plot. 

Scalar  diffraction  theory  predicts  the  existence  of 
on-axis  zeros  in  irradiance  in  the  intervals  predicted  by 
Equation  (14).  This  is  the  first  on-axis  zero  seen  in 
Figure  5  on  the  aperture  side  of  the  focus.  The  theoreti¬ 
cal  irradiance  cross-sections  which  include  this  first 
on-axis  zero  appear  in  Figure  9  for  the  four  cases  under 
study,  and  phase  is  shown  in  Figure  10.  The  two-dimensional 
centerline  plots  were  given  in  Figure  1  in  Section  I.  As 
can  be  clearly  seen,  the  on-axis  zero  disappears  when 
aberrated  with  -0.8  waves  of  spherical  aberration.  Addi¬ 
tionally,  apodisation  removes  the  on-axis  zeros  as  well  as 
reducing  the  amplitude  of  the  modulus  over  the  entire  struc¬ 
ture.  The  apodisation  also  reduces  the  intensified  ringing 
caused  by  the  aberration. 


The  Fresnel  Region .  The  region  centered  on  z/f 
=  0.95  was  studied  using  radial  plots  from  the  optical  axis 
outward  to  normalized  radial  distance  of  50.0.  The  z-axis 
was  taken  from  the  13th  on-axis  zero  from  the  focus  to  the 


.  9.  Cross-sectional  Irradiance  at  the 
First  On-axis  Minimum 
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Fig.  10.  Cross-sectional  Phase  in  Radians 
at  the  First  On-axis  Minimum 
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Fig.  10 — Continued 
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where  n  is  the  number  of  the  on-axis  zero  desired.  In  this 
case  the  values  of  z/f  where  n  =  13  and  n  =  10  define  the 
boundaries  of  the  region.  The  irradiance  plots  of  the  four 
considered  cases  in  this  region  are  in  Figure  11  and  the 
corresponding  phase  plots  are  in  Figure  12.  There  are  some 
significant  predictions  in  these  plots.  The  unapodised 
unaberrated  plot  predicts  an  addition  of  a  fine  wave  struc¬ 
ture  riding  on  the  Fresnel  pattern.  This  fine  structure 
seems  to  intensify  as  the  radius  increases.  The  aberrated 
plot  predicts  a  more  extreme  radial  pattern  with  a  more 
intense  outer  ring.  Additionally,  the  fine  structure  is 
slightly  reduced  in  the  area  outside  the  outer  ring.  The 
fine  structure  has  a  period  corresponding  to  that  given  by 
Young's  double  slit  interference  with  the  slits  taken  at 
the  aperture  edges.  The  apodised  plot  indicates  what 
appears  to  be  a  complete  loss  of  ringing.  On  the  aberrated 
plot,  the  apodisation  again  is  seen  to  substantially  reduce 
the  ringing. 

The  first  minimum  cross-sectional  plane  and  the  12th 
are  studied  experimentally  in  the  next  section  and  provide 
interesting  comparisons  with  the  computer  predictions. 
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IV .  Experimental  Results 
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Experimentation  was  undertaken  to  study  the  effects 
of  apodisation  and  to  test  the  theoretical  predictions. 

The  cross-sectional  planes  studied  were  at  the  first  and 
twelfth  minima  on  the  optical  axis  from  the  focus  toward 
the  aperture.  In  the  unaberrated  and  unapodised  system, 
these  minima  correspond  to  zeros  in  irradiance  and  lie  in 
the  meridional  plots  of  the  previous  section.  The  system 
used  experimentally  consisted  of  a  He-Ne  laser  producing  a 
coherent  beam  with  a  wavelength  of  632.8  nm.  The  aperture 
diameter  was  set  at  1.3335  cm  (0.525  inches)  and  the  system 
was  set  up  to  have  an  f-number  of  12.  An  apodiser,  designed 
and  constructed  by  Mills  (21:178)  and  having  the  character¬ 
istics  given  in  Figure  3  was  used.  The  list  of  experimental 
equipment  is  in  Table  II  and  was  arranged  as  depicted  in 
Figure  13. 

The  system  was  aligned  using  the  techniques  of 
Taylor  and  Thompson  (24:440)  and  Mills  (21:111-115). 

Detailed  explanations  of  alignment  techniques  are  explained 
in  these  sources  and  are  not  covered  here. 

Measurements  of  Aberration 

The  measurement  of  the  amount  of  spherical  aberra¬ 
tion  present  within  the  system  was  accomplished  by  use  of 
a  Smartt  Point  Diffraction  Interferometer  (PDI)  (23).  This 
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TABLE  II 


EXPERIMENTAL  EQUIPMENT 


Equipment 

Model 
or  Type 

Source 

Description 

Laser 

He-Ne 

2mW-linearly  polarized 

Absorbing 

Filter 

945A 

Newport 

Variable  values 

Mirrors 

GM-2 

Newport 

5cm  diameter,  pyrex 

Spatial 

Filter 

LPSF-10  0 

Jodom 

10X  microscope  objective 
10  um  pinhole 

Lens  1 

Aeger s 

80mm  diam.  f  =  495mm 

Lens  2 

Aegers 

54mm  diam.  f  =  580mm 

Iris 

ID-1.5 

Newport 

18  curved  leaves 

Apodiser 

#6 

Mills 

As  in  Figure  3 

PD  I 

Ealing 

Smartt  Point  Diffraction 
Interferometer 

Camera 

4x5 

Graphic 

Poloroid  Type 

Camera 

K1000 

Pen tax 

35mm  used  without  lens 

Film 

57 

Polaroid 
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was  placed  at  the  focus  and  allowed  the  beam  to  pass  through 
the  PDI  while  a  nonaberrated  expanding  wave  front  was  pro¬ 
duced  in  the  focal  plane.  The  pattern  resulting  from  the 
interference  of  the  two  waves  indicated  the  amount  of  aber¬ 
ration  present.  Aberrations  could  be  induced  within  the  sys¬ 
tem  by  reversing  Lens  L2  from  its  depiction  in  Figure  13. 
Since  Ll  and  L2  are  achromatic  doublets,  they  produce  a 
minimum  amount  of  spherical  aberration  when  the  infinite 
conjugate  sides  of  each  lens  face  one  another.  When  L2  is 
reversed  a  maximum  amount  of  negative  spherical  aberration 
is  created.  Zero  aberration  was  measured  when  the  con¬ 
jugates  faced  each  other.  The  f/12  aperture  on  the  PDI 
was  used  and  gives  the  sharpest  contrast  in  fringes  when 
used  with  an  f/12  system.  The  amount  of  aberration  was 
measured  by  comparing  photographs  of  the  image  of  the 
aperture  when  fully  open  and  when  set  at  0.525  inches.  The 
center  of  the  first  dark  fringe  occurs  when  the  interfering 

waves  are  out  of  phase  by  a  half  wavelength.  Since  spheri- 

4 

cal  aberration  increases  radially  as  p  ,  the  aberration  is 

greater  when  the  aperture  iris  is  fully  opened.  When  the 

iris  is  set  at  the  desired  aperture  size,  the  beam  is 

truncated  and  the  amount  of  aberration  is  reduced.  The  two 

photographs  were  then  compared  with  a  simple  plot  of  on 

an  x-y  plot.  The  first  maximum  corresponds  to  one  wave  of 

aberration.  The  truncated  aberration  then  simply  is  taken 
4 

from  the  p  plot  as  m  Figure  14.  The  aberration  was 
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measured  at  -0.78  +  0.05.  The  value  of  -0.8  was  used  in 
the  running  of  all  the  computer  codes  in  which  a  non-zero 
value  for  spherical  aberration  was  needed. 

The  First  On-axis  Minimum 

Comparisons  of  the  four  cases  under  study  at  the 
first  on-axis  minimum  are  in  Figures  15a-d,  along  with  the 
corresponding  computer  predictions  for  the  radial  irradi- 
ance . 

The  alignment  in  this  region  is  extremely  sensitive 
making  the  production  of  a  perfectly  symmetric  diffraction 
pattern  quite  difficult.  This  could  be  caused  by  other 
existing  low  amplitude  aberrations  produced  by  the  test 
lens  L2 .  These  aberrations  were  not  large  enough  to  be 
measured  by  the  PDI  using  the  technique  just  described. 

The  results  of  the  experimental  work  shows  quali¬ 
tatively  good  comparison  with  the  theoretical  predictions. 
The  unaberrated  on-axis  zero  is  lost  due  to  the  loss  of 
edge  diffraction  by  the  aperture  when  the  apodiser  is  in 
place.  Without  the  apodiser,  the  aberrated  system  produces 
a  minimum  at  the  same  axial  point,  though  the  irradiance  is 
greater  than  zero.  In  all  cases,  the  apodisation  is  seen  to 
dramatically  reduce  the  diffraction  ringing  effect.  Addi¬ 
tionally,  there  is  also  a  reduction  in  the  irradiance 
amplitude . 


UNABERRATED  UNAPODISED 


R 


Photograph  and  Calculated  Irradiance 
for  the  First  On-axis  Minimum.  Case  I 
Unaberrated,  Unapodised 


The  Fresnel  Region 

The  12th  on-axis  minimum  was  photographed  for  all 
four  cases.  The  results  appear  in  Figures  16a-d,  along 
with  computer  predictions  for  the  radial  irradiance. 

The  on-axis  minimum  has  an  irradiance  value  of  zero 
for  both  the  aberrated  and  unaberrated  systems.  The  radial 
ringing  has  higher  peaks  and  lower  valleys  for  the  aber¬ 
rated  system  which  are  seen  as  brighter  and  darker  regions 
in  the  photograph. 

The  radial  fine  structure  is  clearly  seen  on  both 
unapodised  photographs.  They  show  an  angular  dependence 
on  this  structure  due  to  a  non-circular  aperture.  The 
aperture  used  was  an  iris  with  18  curved  leaf  edges  arranged 
in  a  circular  pattern.  At  optical  wavelengths,  this  circu¬ 
lar  approximation  is  easily  detected  by  the  angular  depen¬ 
dence  of  the  fine  structure.  Additionally,  the  region  near 
the  optical  axis  shows  angular  dependence.  This  is  also 
due  to  the  non-circular  aperture. 


\V 
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ABERRATED  UNAPODISED 
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Fig.  16b.  Photograph  and  Calculated  Irradiance 

for  the  Twelfth  On-axis  Minimum,  Case  II 
Aberrated,  Unapodised 


UN  ABERRATED  APODISED 


so 


Photograph  and  Calculated  Irradiance 
for  the  Twelfth  On-axis  Minimum,  Case  III: 
Unaberrated,  Apodised 
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ABERRATED  APODISED 


Fig.  I6d.  Photograph  and  Calculated  Irradian' 
for  the  Twelfth  On-axis  Minimum,  C 
Aberrated,  Apodised 
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V.  Conclusions 


The  theoretical  predictions  and  experimental  evi¬ 
dence  clearly  demonstrate  that  diffraction  effects  are  sub¬ 
stantially  reduced  when  an  aperture  is  apodised.  The  ring¬ 
ing  phenomenon  is  completely  attenuated  in  the  unaberrated 
case  and  substantially  reduced  in  the  aberrated  case.  Also 
the  irradiance  is  also  reduced  in  all  cases  of  apodisation. 

Additional  conclusions  of  this  research  follow: 

1)  Scalar  diffraction  theory  provides  enough  accu¬ 
racy  necessary  to  study  this  phenomenon, 

2)  The  paraxial  and  Fresnel  approximations  must  be 
carefully  considered  in  the  region  of  study  so 
as  to  not  limit  the  accuracy  of  the  theory, 

3)  At  optical  frequencies  and  aperture  sizes  on 
the  order  of  a  centimeter,  the  amount  of  data 
required  to  model  near-field  irradiances  becomes 
extremely  large, 

4)  When  spherical  aberration  is  present,  the  normal 
Fresnel  diffraction  patterns  are  more  intense  in 
the  peaks  and  less  intense  in  the  valleys, 

5)  The  predicted  fine  structure  in  the  Fresnel 
region  is  real  and  was  observed  experimentally; 
additionally,  as  spherical  aberration  increases, 
the  fine  structure  decreases.  This  structure 


has  the  same  amplitude  and  frequency  as  that 
matching  Young's  double  slit  interference  with 
slits  at  the  aperture  edges. 

Areas  for  Further  Study 

The  ability  to  model  the  near-field  exists  with  the 
use  of  the  codes  in  this  study.  The  cross-sectional  plane 
can  be  accurately  predicted  by  the  code  N256.FOR  which  can 
handle  all  third-order  aberrations  and  all  Gaussian 
apodisers.  The  meridional  plane  can  accurately  be  predicted 
by  the  use  of  IRAD.FOR,  which  can  also  handle  non-circularly 
symmetric  aberrations.  The  only  limitations  on  the  use  of 
these  codes  is  the  computer  hardware  and  the  regions  where 
the  paraxial  and  Fresnel  approximations  maintain  good  accu¬ 
racy.  With  optical  frequencies  and  large  f-numbers  this  is 
a  significant  portion  of  the  Fresnel  region. 

This  study  should  be  continued  in  depth  with  the 
study  of  the  other  third-order  aberrations,  though  this 
would  be  more  difficult  experimentally.  Additionally,  this 
study  should  be  repeated  with  lasers  of  high  energy. 

Also,  a  more  complete  study  of  phase  should  be  made. 
The  treatment  of  phase  incorporated  in  this  study  is  only  a 
cursory  treatment  and  it  should  be  investigated  more 
thoroughly . 


Appendix:  Computer  Codes 


Several  computer  codes  were  written  for  this  thesis 
of  which  three  are  included  in  the  appendix.  These  codes 
include  the  most  important  programming  techniques  employed. 
The  plots  in  the  body  of  the  thesis  were  drawn  with  the 
software  package  DISSPLA  (6)  and  printed  on  a  laser  plotter. 


PROGRAM  N256 


C 

C  THIS  program  calculates  the  phase,  modulus  and 

C  INTENSITY  OF  THE  TWO-DIMENSIONAL  IMAGE  OF  AN 
C  ON-AXIS  POINT  SOURCE  THROUGH  AN  OPTICAL 
r.  SYSTEM  WITH  ARBITRARY  ABERRATIONS  AND 
C  APOD  I  SAT  I  ON  INCLUDING  THE  NEAR  FIELD 
C 

c 

I NTEGER  I Wf  < 2000 )  .  I A 1 ,  I A2 . N 1 , N2 . N3 , I JOB 
INTEGER  S2.S 
REAL  RW<  <2000 

REAL  c  '  CCS.  25c;;  .  M  (256.  25  c- ;  .  E  (256)  .  I  NT  <256, 250 
COMPLEX  A  <256. 256?  . CWt  <256? , B  <256. 256' , C  <256? 
COMPLEX  D (256? 

WRITE  2  NC'UT  S:’ 

READ  <*.#•  E 

WRITE  *.*?  ’INPUT  RADIUS  MULTIPLIER:’ 

READ  <*,*)  RM 

WRITE  <*.*)  ’NORMALIZE?  <1.0  =  YES,  0.0  =  NO):’ 
READ  <  * . * )  FLAG 
IF  (FLAG  . NE.  1.0)  THEN 

flag  =  o.o 

W,R  2  TE  (  *  ,  '4  ’  -Jn r< or  .ns.l  ;  r  sc  ' 

ENE  Ir 


OPEN (UNIT=22. STATUS=’ NEW’ , FILE=’ PH. DAT' ) 
OPEN  <UNIT  =  23. STATUS=’ NEW’  , P I LE= ’ MOD . DAT ’ > 
OPEN (UNI T=24. STATUS^’ NEW’  , F I LE= ’  I NT.  DAT’  ) 
OPEN (UNI T  =  25. STATU£=’ NEW' . F I LE= ’ BB . DAT ’  ) 

S2  =5/2 

C 

C  c  EAD  PARAMETER  VALUES  FROM  MIC  FILE 


C  READ <5. 20) N 

N  =  256 

C  READ (5. 40) A20 

A 20  =  0. 

WRITE (*.*)’ A21 : ’ 
READ ( * . *) A21 
C  A21  =  0.01 

C  READ <5. 40) A22 

A22  =  0. 

:  REAL- :  5. 40?  AC  1 
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n  n  o 


A3 1  =  O. 

C  READ <5. 40) A32 

A32  =  0. 

WRITE  (*.*>  'INPUT  A42: ' 

READ  <  *. *  > A42 

WRITE  <*.  *>  ' INPUT  GG: ' 

READ ( t . *)GG 

WF.  I TE  <  # , #  >  'INPUT  z.'-f:’ 

READ  '  »  .  *  >  FRAC 

C  WRITE (5. *  >  'INPUT  APERTURE  FRESNEL  NO.:  ' 

2  READ  (5.#)  FN 

FN  =  439.021 
PI  =  4 . *  AT AN (1.0) 


2  20  FnRMAT';J: 

40  FORMAT ( 1FE1 1 . 4) 

41  FORMAT < 14. 3X. 14. 3X , 1PE1 1 . 4, 3X . 1PE1 1 . 4) 

42  FORMAT <14. 3X.  I  4 . 3 X .  1  PE  1 1 . 4 ) 

a-  5rZ'PMA"r  ■  3  '2X .  I  -2  *  .  Z  '3X .  PS.  4 )  > 

FOPMA,T  14. 3  >. .  1  P£  :  ;  .  4 . 3  >  .  1  PE  1 1 . 4  > 

r 

C  N=NUMBER  OF  SAMPLE  F 01 NTS  IN  ONE  DIMENSION 

C  A20=4w  deo  ASTIGMATISM 

L  A21=F0CAL  SHIFT 

2  A 2 2=0  deg  ASTIGMATISM 

C  A3 1 =X  COMA 

2  A3Z-Y  COMA 

C  A42'=SF'HER'ICAl 

2  GG-WIDTH  0~  GAUSSIAN 

C  FF AC=F'ERCENT AGE  OF  FOCAL  LENGTH  IN  NEAR  FIELD 

2  FL=FOCAL  LENGTH  IN  APERTURE  DIAMETERS 


C 

r 

V- 

C  CIRCULAR  APERTURE  WITH  VARIABLE  TRANSMITTANCE 
2  AND  ABERRATIONS 
C 


PI  =  17.0#RM 
P'2  =  1E.0#RM 


DO  220  I  =  1  .  N 

DO  200  J  = 1 . N 

RAD=SORT ( < I  -  (R2+0. 5) ) **2+ ( J- <R2+0. 5) ) **2) 
X= ( I- (R2+0.5) ) /R1 
Y  = (J- (R2+0.5) ) /R 1 

CALCULATE  CONTRIBUTIONS  OF  THIRD  ORDER 
ABERRATIONS 

USING  ZERNIKE  MONOMIAL  REPRESENTATION 
ARG20=A20*Z.  *>:*  V 
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p 


p 


5 


ARG2 1 =A2 1 t (-1+2. *Y*Y+2. *X*X ) /2. 

AF:G22=A22*  <Y*Y-X*X ) 

ARG3 1 =A3 1  * (-2.  *X  +  3. *X*Y*Y+3. *X**3) /3. 

ARG32=A32* (-2. *Y+3. *Y**3+3. *X*X*Y> /3. 

ARG42A  =  A42*  <  1.  -  6. *Y*Y  -  6.*X*X  +  6.*Y**4  ) 

ARG42  =  ARG42A  +  (  12. *X*X*Y*Y  +  6. *X**4  >  *A42 
ARG42  =  ARG42/6. 

ARG 1  =  ' ARG20+ARG2 1 +ARG22+ARG3 1 + ARG32+ARG42 >  * 2 . *P I 
C  CALCULATE  CONTRIBUTIONS  OF  VARIOUS  AF'ODISERS 

C  GAUSSIAN  ARODISER 

C  TT  =EXF'  <-GG*RAD*RAD/Rl  /R1 ) 

C  SUPERGAUSS I AN  ARODISER 

F=  1  . 

R  =  RAD  /  R 1 

TT  =  E  >  F  i-GG  *  F:U2) 

C  CALCULATE  ARG2 .  THE  NEAR-FIELD  CONTRIBUTION 

ARG2  =  FI  »  FN  *  (1  -  FRAO  /  FRAC  *  R**2 
I F  (RAD . LT . R3+0. 5) GO  TO  60 
1“  (PAL  .  c~.  F.1+0.  5'  GO  TC  100 
c,c-=-c.AD*F' 1  +  0.5 
GO  TO  80 

t  :  a • I . J> =CMPLX (TTtCOS (ARG1+ARG2) , TT*S IN < ARG 1+ARG2 > > 

GO  TC  120 

60  A (I. J)=CMPLX (SC#TT#C0S (ARG1+ARG2) , SC*TT*SIN ( ARG1+ARG2) ) 

GO  TO  120 

i  f.'-.1  A  1  I  .  3  )  —  ( O  •  .  V .  ) 

: 2 :  CONT I NUE 

2..  CONTINUE 

22!  CONTINUE 
0  DO  230  1=1.  s 

Z  DO  225  J=1.S 

0  M ( I . J  > =CABS  < A ( I . J ) > 

C  WRITE <24. 42) I .J. MCI. J> 

C  215  CONTINUE 

C  23:  CONTINUE 

C 

C  PERFORM  FAST  FOURIER  TRANSFORM  USING  1MSL 


240  I A  1  =r  J 
IA2-N 

f  ;  -  - 

N3=  1 
I J0B=-1 

CALL  FFT3D  (A.  IA1 , IA2, N1 . N2, N3, I  JOB, IWK, RWK, CWK) 
C  DO  242  J=i , N 

C  1  =  1 

C  WF.  I TE  ( 22 .41)  I ,  J ,  A  ( I ,  J  ) 

C242  CONTINUE 
C 

C  REDUCE  SIZE  OF  ARRAY 


DO  2 SO  1=1. S2 

DO  245  J=1 . S2 

B ( I  - J)=A(1, J) 

B ( I . S2+ J  >  =A  (  I . N-S2+J ) 

8 ( S2+  I . J  >  =  A (N-S2+ I . J) 

B (E2  + I . S2+ J ) = A <  N-S2  + I . N-S2+ J ) 
245  CONTINUE 

250  CONTINUE 

C  DC  252  1=1.5 

C  DO  251  J=1.S 

C  WF: I TE  (23.41)  I .  J  ,  B  ( I .  J  > 

C25 1  CONTINUE 

C252  CONTINUE 

C 

C  :.,:t  2  ERG  FREQUENCY  TO  ARRAY  CENTER 


DO  260  1=1.5 

DC  255  J= 1 . S2 

A  <  I .  J  >  =5  '  I  •  S2+  1 —  J 
a  <  i . S2+J  > -b: : . s*:-J) 

CONTINUE 
CONTINUE 
DO  270  1=1,52 
DO  265  J=1,S 

B  ( I  , J) =A (S2+1-I . J ) 
B(S2+I. J)=A<S+1-I,  J) 
CONTINUE 
CONTINUE 
DC  2-4  J  =  1 .  £• 

DC  272  1=1.5 

TYRE  41, I  .  J  ,  B  C I ,  J  > 
CONTINUE 
CONTINUE 


C  FIND  PHASE  OF  ARRAY 

f 

U 

243  DO  360  I = 1 . S 

DO  340  J=1.S 

AA=A1MAG (B  <  I . J)  ) 

BP  =  RE AL <B  < I . J  )  > 

1  F  •  BE  .  EC; .  •:> .  OR 000000 )  BB=0 .  00000 1 
r  . :  .  J  > =A7AN2 <AA. BB; 
r  WRITE (25,  42) I , J . P ( I , J ) 

34:  CONTINUE 

36''.'  CONTINUE 

c 

C  FIND  MODULUS  AND  INTENSITY  OF  ARRAY 
C 

DO  400  I = 1 . S 

DC  380  J  =  1  .  S 

M(I.J)  =  CAPS  (  B  (  I , J )  )  /  FRAC 
I  NT  l I  . J ) =M ( I . J)  **2. 


u  U  u 


380  CONTINUE 

PRINT  *. ' INT(S2. S2)  =* , INT  <S2, S2) 
CENT  I  NT  =  INT (S2. S2) 

CENTM  =  M(S2.S2> 

400  CONTINUE 

IF  (FLAG  . EQ.  0.0)  GOTO  401 
DO  21  1=1 . S 

DO  22  J= 1 . S 

M(  I . J) -M ( I , J) /CENTM 
INT  <I.J)=INT ( I , J ) /CENT I NT 
22  CONTINUE 

3 '  CONTINUE 
401  CONTINUE 

PC'  404  J  =  1  .  S 

DO  402  J  =  1 . S 

WRITE  <23. 42)  I.J.M(IJ) 
WRITE ( 24. 42 ) I . J . INT ( I , J ) 
CONTINUE 
CONTINUE 

c  WRITE  CENTRAL  SLICE  OF  INTENSITY 
DO  415  I = 1 . S 

C  WF:ITE<23.44)  I.MU.S2) 

WRITE (25. 44) I, INT(I.S2) 

413  CONTINUE 

DL  =  3. *PI/2. 

0!  =  3.  *  P  I  /  4  . 


REMOVE  TILT  FROM  EACH  ROW 

DO  420  1=1 . S 
DO  410  J=2.S 

P< I ,  J) =P(I , J) - (J-l> *.4295 
WRITE <23, 42) I , J , P  < I , J ) 

410  CONTINUE 

4  CO  CONTINUE 


DO  500  I=1,S 
DEL=0. 

I CNT =0 

DO  490  J=1.S 
I CNT  = I CNT  + 1 
PCI. J)*P(I, J)+DEL 
F 1 =P ( I , J ) 

G 1 =P ( I . J  + 1 ) +DEL 
H 1 = ABS (01  — F 1 ) 

WRITE (22. 43' I. J. I CNT . DEL , F 1 . G 1 .HI 
IF (HI . LE. Q1 ) GO  TO  490 
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U  U  U  U  U  i  L  j  U  <_J  LJ  U 


IF( (Hl.GE.DL) .AND. (Gl.GT.Fl) ) GO  TO  450 
IF( (Hl.GE.DL) .AND. (G1.LT.F1) ) GO  TO  460 
IF ( ICNT . LT. 5) 60  TO  411 
I CNT  =  4 

C  DO  LINEAR  FIT  TO  LAST  ICNT  POINTS (B1=SL0PE> 

41 :  IF < ICNT.EO. 1 >G0  TO  431 

SX Y=0 . 

sx=o. 

SY— 0. 

SX2=0. 

DO  421  K=1 . ICNT 
X  =  I CNT  + 1  — K 
SXY=SXY+X*P ( I . J+l-K) 

£>'rC.y  +  | 

SY=SY+P ( I . J+l-K) 

SX2=SX2+K**2. 

4  2 :  CONT I NUE 

B 1  =  f  SXV-SX  *SY  / 1  CNT  ■  /  <SXZ-SX*SX  ‘ICNT) 

GO  TO  42- 

471  Bl=PU.J+2)-P(I,J+l) 

42-  1F( (Bl.LE.O. ) .AND. (Gl.GT.Fl) ) GO  TO  430 

IF( (Bl.GT.O. ) .AND. (G1.LT.F1) ) GO  TO  440 
I CNT =0 
GO  TO  490 

47  0  I NCT =0 

GO  TO  450 

44;  I CNT  =0 

GO  TO  460 

42:  DEL=DEL-2. *PI 

GO  TO  490 

4-:  DEL  =  DEL  +  2.  *F'I 

4°:  CONTINUE 

507  CONTINUE 

C  DO  850  I=1.S 

DO  849  J=1.S 

WRITE (24, 42) I , J , P ( I , J ) 

64=  CONTINUE 

S3:  CONTINUE 

:  *•  ~  Z0T  RM4.SE  ON  entire  ROWE  RELATIVE  ”?  COW  I=* 

DO  1 1 60  I = 1 . S 
J=S2 

WRITE (25, 42) I . J , P  < I . J ) 

11=7’  CONTINUE 

REMOVE  TILT  BETWEEN  ROWS 
DC  1200  1=2. S 

DO  1150  J=1.S 

P ( I . J ) =P  ( I , J ) -  ( I - 1 >  * . 4295 


C  J  =  S2 

C  WRITE  (22.  42)  I.J.F'(I.J) 

Cl 210  CONTINUE 

C 

C  HANDLE  PHASE  JUMPS  BETWEEN  EACH  ROW 
C 

DEL=0 . 

ICNT=0 

DO  1500  I = 1 . S 
I CNT = I CNT  + 1 
DO  1510  J=1.S 
P ( I , J ) =P  C I . J ) +DEL 
J51 :  CGMTINUE 

1SC'' 

F 1 =P ( I . J ) 

G1=P ( 1  +  1 . J  > +  DEL 

M 1  -ABE  ■'  G 1  —  F  1  ■ 

;  NF  I  TE  '24.471  I  .  J  .  I  CNT  .  DE_  .  r  1  .  G 1  .  H 1 

IF (HI. LE. 01) GO  TO  1500 

IF ( (Hl.GE.DL) .AND. (6i.GT.Fl) ) GO  TO  1450 
IF< (Hl.GE.DL) .AND. (Gl.LT.Fl) ) GO  TO  1460 
IF ( ICNT. LT. 4) GO  TO  1411 
I CNT =2 

DC  LINEAR  FIT  TO  LAST  ICNT  POINTS 
1411  IF ( ICNT.EO. 1 >60  TO  1426 
S.<  v=0 . 

S  v  =  0 . 

SY=0. 

SX2=0. 

DO  1421  K= 1,1 CNT 
X  =  I CNT  + 1  — K 
SXY=SXY+X*F ( I + 1 -K . J > 

SX*SX+K 

SY=SY+F  ( I  +  l-K. J  > 

SX2=SX2+K**2. 

142.:  CONTINUE 

&1=<SXY-SX*SY/ICNT> / (SX2-SX*SX/ICNT> 

GO  TO  1427 

;  .  El  =P  '  T  MZ .  J  '  -P '  I + 1 . 2  • 

if  .  •:  E 1  •  t.  E .  0 .  )  .  AND .  (E..GT.F1'  *60  T0  1  4  70 

IF ( CB1 . GT. 0. ) . AND.  <G1 . LT. FI ) ) GO  TO  1440 
I CNT =0 
GO  TO  1500 
1470  I CNT =0 

GO  TO  1450 
1440  I CNT  =  0 

GO  TO  1460 

1450  DEL=DEL-2. *PI 

GO  TO  1500 

14*'  DEL=DEL+2.  *F-I 

1500  CONTINUE 

DO  2470  J  =  1.S 
I=S2 

C  WR I TE ( 25 .42) I.J.P(I.J) 

24-'0  CONTINUE 

C  CENTEF;  PHASE  AROUND  ZERO  VALUE 


% 


♦:V 


’VJ 


4i 


/»«| 


m 


f|*i| 

5jl; 

r^i«. 


y>l 

& 


$ 1 


fid 


i 


,vj 


DO  250 1  1=1 . S 
DO  2491  J=1 . S 

F'<I.  J)=P  (I,  J)-9. 

WRITE (22.42) I , J,R  U . J> 
CONTINUE 
CONTINUE 


C  WRITE  CENTRAL  SLICES  OR  MODULUS  AND  PHASE 


DC  2550  J=1.S 

WRITE (22. 44) J. P ( I . J) 
WRITE (25. 44) J. M<I. J) 
CONTINUE 
CLOSE ( UN I T =22 ) 

CLOSE (UNIT=2?) 

CLOSE (UNIT=24) 

CLOSE  <UNIT=25) 

STCF 

END 
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PROGRAM  MERID.FOR 

C  MERID  CALCULATES  THE  NEARFIELD  IRRADI ANCE  IN  THE 

C  MERIDI ANAL  PLANE.  IT  A3SUMS  A  LENS  AT  THE  APERTURE  WITH  4 /# 

SET  AS  INPUT.  ADDITIONALLY  ’ MERID. FOR'  TAKES  INTO  ACCOUNT 
THE  CIRCULARLY  SUMMETRIC  ABERRATIONS  OF  DEFOCUS  AND  3RD  ORDE 

C  SPHERICAL  ABERRRATION.  THE  MODULUS.  PHASE.  AND  THE  OPTICAL  A 

C  IRRADI ANCE  IS  ALSO  CALCULATED. 

PARAMETER  (PI  =  3.141592654) 

REAL  IRR (256. 256) . M (256. 256) . PH (256. 256)  , MMESJO. YI (256) 

REAL  Yp (256) 

INTEGER  I.J.K.KK 

CHARACTER* 14  F I LNAM 1 , F I LNAM2 , F I LNAM3 

2  data 

flag  =  l.o 

WRITE  (*.*)  ’OPTICAL  AXIS  FILE'"'  ("1.0"  FOR  YES.  "0.0"  FOR  NO 
READ  <*.*)  FLAG 1 

WRITE  <*.*)  ' INPUT  IRR  FILENAME:’ 

READ  (*.2)  FILNAM1 

WRITE  (*.*)  ’  MOD  FILENAME:’ 

READ  <  * . 2 '  FILNAM2 

WRITE  (*.*)’  PH  c I LNAME .  ("0"  FOR  NULL':’ 

READ  (*.2)  FILNAM3 
IF  (FILNAM3  .EG).  ’O’)  FLAG  =  0.0 
2  FORMAT <  A 1 4) 

1  FORMAT (F10. 4) 

WRITE  <*. *)  ’ INPUT  A21 : ’ 

READ  (*,*)  A21 
C  A2 1  =  0 . 0 

WRITE  (*,*)  ’INPUT  AA2:' 

READ  <*,*)  A42 
C  A42  =  0.0 

WRITE (*.*)  ’ INPUT  GG: ’ 

READ  '*.*'■  Gr- 
56  =  0.0 
FN  =  439.021 

WRITE <*.*)  ’INPUT  START  C: ’ 

READ <  *  .  * )  Cl 

WRITE <*.*)  ’INPUT  END  C:  ’ 

READ  (*.*>  CEND 

WF.’ITE  (  * .  * )  ’INPUT  START  X:’ 

READ  (*.*>  XO 

WRITE (*.*>  ’INPUT  END  X:’ 

READ  <*,*)  XEND 

WF I TE ( * . * )  ’ ARRAY  ELEMENTS  IN  Z  DIRECTION:' 

READ  ( * . * )  LASTC 
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WRITE  <*.*>  ' ARRAY  ELEMENTS  IN  R  DIRECTION:’ 

READ  <*.*)  LASTR 
R 1  =  XO 

WRITE  (*.*)  ’ INPUT  KK.  (must  be  an  odd  integer):  ’ 

FORMAT (14) 

READ (  *  .  7 )  KK 
PRINT  *. 'Ll  ='  .  KK 

DC  =  (CEND  -  Cl) /(FLOAT (LASTC)- 1.0) 

DF.  =  (XEND-XO)  /  (FLOAT  (LASTR  >-1.0) 

PRINT  *. ’ DR  =  ’ , DR 

4  FORMAT ( 1 X. 9X. ’ J=’ . 13) 

3  FORMAT (IX. ' I = ’ , 13) 

F  =  R 1 
C  =  Cl 

DRHO  =  1.0/ (FLOAT (KK) -1.0) 

PR  I  NT  # .  ' DRHD= ' . DRHO .  ’ LASTR= ’  , LASTR 
DC  10  1=1. LASTC 
WFITE  (*.7)  3 

DO  20  J=1 . LASTR 
C  =  I*DC  +  Cl 
R  =  (J-l ) *DR  +  R 1 
WRITE  (*.4)  J 
RHO  =  0 . 0 
DO  30  K= 1 . KK 

*  CFEATjrJS  THE  INTEGRAND  * 

*  Aberration  terms  using  Ze’-nike  F'ol  vnomi  al  s  * 

ARG21  corresponds  to  de-focus 

ARG42  corresponds  to  spherical  aberration 

ARG21  =  A21  *  <2.0  *  RH0**2  -  1.0) /2.0 

ARG42  =  A42  *  ( 6. .  0*  < RHO* *4  -  RH0**2)  +  1.0)  /6.0 

ARG1  =  (ARG21  ♦  ARG42) *2. 0*PI 

*  Aoodisation  Term  * 

TT  =  E  /  F  ’  -GG*PHC*  *  2 ' 

*  Fresnel  Region  Term  * 

ARG2  =  PI  *  FN  *  (1-C)/C  *  RHO* *2 

REA  =  TT*C0S (ARG1+ARG2) 

AIM  =  TT*SIN (ARG1+ARG2) 

ARG  =  PI  *  R  *  RHO/C 

YR  < K )  =  REA  *  RHO  *  MMBSJO (ARG. IER) 

IP  (IEF  .  NE.  0)  PRINT  *  .  '  IEC  =  '  .  IEF; 

Y I  < K )  =  AIM  *  RHO  *  MMBSJO (ARG.  IER) 

IF  (IER  .NE.  0)  PRINT  * . ' I  IER=' , IER 


RHO  =  RHO  +  DRHD 
7<;>  CONTINUE 

ODDR  =  0.0 

ODDI  =  0.0 

DO  40  K=3.Kk— 2.2 

ODDR  =  ODDR  +  2.0tYF:(K) 

ODDI  =  ODDI  +  2 .  O  #  Y I ( K ) 

4<;  CONTINUE 

EVENR  =  0.0 
EVEN I  =  0.0 
DO  50  k'=2.  r-'k.-l  .  2 

EVENR  =  EVENR  +  4.0»YR<k) 

EVENI  =  EVENI  +  4.0*YI(K) 

50  CONTINUE 

ElN'r  =  i  DR  HO  3.  O'  *  '  vc‘  '  1  >  *  EVENK  +  ODEc  +  vp  '■*'#  >  > 

All  NT  -  (DRH0/3.0)  *  <  Y I  ( 1  >  +  EVENI  +  ODDI  +  YI  <KK)  ) 

M'l.J)  =  2.  O-'C  *  SORT  (RINT#*2  *  AIINT**2' 

I RR  < 1.3)  =  M  <  I .  J )  *  *  2 
PH(I.J)  =  AT  AN2 (All NT . R I NT ) 

20  CONTINUE 

10  CONTINUE 

OPEN (UNI T=22 . ST  ATUS= ’ NEW 7  , FILE=FILNAM1 ) 

IF  (FLAG!  .  EO .  1.0)  OPEN  (UNIT =20.  STATUS='  NEW’  ,FILE=’A>'IS.  DAT '  ) 

OPEN ( UN I T  =  24 , ST ATUS= ?  NEW ■  , F I LE=F I LNAM2 ) 

OPEN  ( UN I T  =25 . ST ATUS= '  NEW’  , F I LE=F I LNAM3 ) 

44  kqkmaT  <13.  3X  .  13.  3X  ,  IF'EI  1 . 4) 

DC  60  I=1,LASTC 

DC  70  J=1 . LASTR 
V  =  IRR(I.J) 

WR I TE  < 22 .  44 )  I.J.V 
WRITE <24. ^4)  I.J.M(I.J) 

IF  •  FL  AC?  .EC.  1.0‘  WRITE  .'25.  44'  I .  J  .  PH  <  J  .  .1 ) 

CONTINUE 
6>0  CONTINUE 

45  FORMAT ( 13. 3X. 1PE1 1 . 4) 

IF  <FLAG1  .EO.  1.0)  THEN 
DO  SO  1=1. LASTC 

V  =  IRR ( I . LASTR/2) 

WRITE <23. 45)  I,V 
80  CONTINUE 
END  IF 

CLOSE  <  UN  I T  =  22 ) 

IF  ( FLAG 1  .EO.  1.0)  CLOSE <UNIT=23) 

CLOSE  <UNIT  =  24) 
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CLOSE  <UNIT=25> 


END 
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PROGRAM  I RAD 

REAL  IRR (2000) . DDL  IN, FL . LAMBDA , M,  K 

INTEGER  I , L 

EXTERNAL  F.AY.BY 

COMMON  GG. K, PI . L , A42, P. Q 

CHARACTER  F I LNAM* 1  A 


F I  =  4 . 0*ATAN ( 1 . O) 


I  no -it 

Parameter  Values 

WF  ITE 

(  *  .  *  ' 

'FILENAME: ' 

FORmAT  'A14> 

READ 

<  * . :  3  > 

p  ilnam 

WR  I  TE 

( * .  * ) 

'X  START:' 

READ 

<* .  *  > 

xo 

WRITE 

<*.  *) 

' X  END:  ' 

R  F"  4-  D 

TEND 

WF  :~e 

'AFRRU  eleme 

READ 

(*.*:>  N 

write 

<  *  .  #  ) 

-  r- .  • 

w  ■ 

READ 

(*,  *> 

C 

LAMBDA  =6.2 

39E-05 

FL  - 

1 6 . 002 

A  =  0.66675 
Z  =  C  #  CL 


WT  7 

T  27 

(*.  * : 

'  A42:  ' 

REA 

D 

it.*) 

A42 

WRI 

TE 

it.*) 

'GG:' 

REA 

r, 

d.*) 

Gb 

DX 

= 

<  XEND 

-  XO) / (FLOAT (N)- 1.0) 

X  - 

XO 

AP 

= 

0 . 0 

BR 

= 

1 . 0 

WF  I 

TE 

(*.*)’ 

AERR:'  1 

PER 

r- 

<  t .  t ) 

AEF.R 

t  - 

- 

.  '*r  I 

'-"re  a 

!•  f  : 

= 

c  ;;«r 

>  i  t  2  ■*  Z  <  X  2  ‘ 

PO 

= 

E*A**2  *  (RRO  -  FL)  /  (2.0  *  RRO  *  FL) 

00  =  K*A*y/RkO 
RP  =  PRO 
D  =  00 
P  =  PO 

4 T  FORMAT < 14. 3X.  1PE1 1 . 4> 

OPEN  CUNrT=l . STATUS—' NEW’  . F I LE=F I LNAM > 


L=  1 

RE  =  DBLIN(F.AR.BR.AY.BY.AERR.ERROR. IER) 

PRINT  *.'L=’.L.'  RE=  ' . RE . '  ERR=  ’  .  ERROR 
IF  ( IER  .  NE.  0  .AND.  65  .AND.  66)  CALL  CHECK < ERROR 
*  IER.I.X) 

L  =  2 

A I  =  DBLIN <F. AR. BR, AY, BY. AERR. ERROR. IER) 

IF  ( I ER  .NE.  0  .AND.  65  .AND.  66)  CALL  CHECK  (ERROR 
4  .IER.I.X) 

PRINT  *.’L=’,L,’  A I  = ' ,  A I , *  EPR=’. ERROR 

m  =  SORT  (RE*# 2  +  AI#*2)  /  (2. 0*PI #FL*RR)  /  1.9526: 

IRR(I)  =  M**2 

WRITE  (1.R2.)  I  .  I  PR  ( I ) 

PRINT  #.’I  =  ?,I.'  IRR='.IRR(I> 

PR  I  NT  * .  ’  ' 

y  —  v  "-t  y 

c  f  =  c.ne  -  •;  >  t  -  7**2> 

p  =  K*A**2* (RR-FL) / (2. 0*RR*FL> 

C  =  K*A*X/RF: 

CONTINUE 
CLOSE (UNI T=1 ) 

END 

***************  Subroutines  ************************ 

F:EA_  F'JN'C’ION  F  (R.  V) 

REAL  LAMBDA. K, PI 
COMMON  GG.K.FI.L.A42.P.Q 

IF  (L  .EO.  1)  THEN 

F  =  EXP (-GG* CR*A) **2>  *  COS(-  P*R**2  -  Q*R*COS(Y>  + 

4  h42»  (R**4  -  R*#  2  *  1 . 0/6. 0) #2*PI )  *  R 

ELSE  IF  (L  . EO.  2)  THEN 

F  =  EXP !-G6* (R*A> **2)  *  SIN<  -P#R**2  -  Q*R*3IN(YJ  + 
5  A42*(R*#4  -  R**2  +  1.0/6.0) *2*PI>  *  R 

END  IF 

::  r  (2-r.i 

LMI. 


SUBROUTINE  CHECK  (ERROR. IER.I.X) 
REAL  ERROR. X 
INTEGER  IER. I 

IF  (IER  .GT.  0)  THEN 

PRINT  *. ’ IER=5 , IER. ’  I=’.I,’  X=7 , X 

PRINT  *.'  ERROR  = ’ . ERROR 


RETURN 

END 

REAL  FUNCTION  AY(R.' 

REAL  R 

AY  =  0.0 

RETURN 

END 

REAL  FUNCTION  BY (R  > 
REAL  P 

BY  =  2.0  *  4. 0* AT AN • 1 

RETURN 

END 
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The  effects  of  Gaussian  apodisation  with  -0.8  waves 
of  spherical  aberration  on  a  coherent  optical,  system  was 
examined.  The  area  of  study  was  in  the  near^field  concen¬ 
trating  on  the  areas  near  the  focus  in  the  Fresnel  region 
of  an  f/12  system  using  optical  wavelengths.  Computer  pre¬ 
dictions  were  made  for  four  cases:  unapodised-unaberrated, 
apodised-unaberrated,  unapodised-aberrated,  and  apodised- 
aberrated.  Predictions  were  made  for  cross}- sectional  planes 
perpendicular  to  the  optical  axis  using  Fourier  optics. 
Meridional  plane  predictions  were  produced  using  a  numerical 
integration  method  of  determining  the  Kirchoff  integral. 
Additionally,  experimental  data  are  given  to  compare  with 
the  predictions.  It  is  shown  that  the  experimental  data 
matches  the  computer  predictions  and  that  apodisation  is  an 
effective  method  for  controlling  the  ringing  due  to  edge 
effects  and  spherical  aberrations.  Additionally,  fine 
structure  corresponding  to  Young's  double  slit  interference 
is  observed  in  unapodised  cases. 
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